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We investigate the detection of exotic massive strongly interacting hadrons (uhecrons) in ultra 
high energy cosmic ray telescopes. The conclusion is that experiments such as the Pierre Auger 
Observatory have the potential to detect these particles. It is shown that uhecron showers have 
clear distinctive features when compared to proton and nuclear showers. The simulation of uhecron 
air showers, and its detection and reconstruction by fluorescence telescopes is described. We deter- 
mine basic cuts in observables that will separate uhecrons from the cosmic ray bulk, assuming this 
is composed by protons. If these are composed by heavier nucleus the separation will be much im- 
proved. We also discuss photon induced showers. The complementarity between uhecron detection 
in accelerator experiments is discussed. 



I. INTRODUCTION 

Ultra high energy cosmic ray (UHECR) observatories 
investigate the high energy end of the cosmic ray spec- 
trum (above ~ 10 19 eV). Their results [31 H] are consis- 
tent with the presence of the Greisen [T] and Zatsepin 
and Kuzmin [5] (GZK) feature. 

GZK showed that nucleons propagating through the 
Cosmic Microwave Background Radiation (CMB) will 
have their energy degraded. The main energy loss mech- 
anism for cosmic rays above ~ 5 x 10 19 eV is pion photo- 
production. In order to reach the Earth, nucleons have to 
be produced relatively near us, at a maximum distance 
of about 100 Mpc. As a consequence, the cosmic ray en- 
ergy spectrum should fall steeply around ~5x 10 19 eV. 
This feature is known as the GZK cutoff. Since there 
are events [6] detected beyond this cutoff, their ori- 
gin, composition and sources became a puzzle and the 
existence of the GZK cutoff was questioned. 

Here we investigate the possibility of detecting exotic 
massive and strongly interacting hadrons (uhecrons) in 
the Pierre Auger Observatory [7J. Uhecrons were first 
proposed [3] as a solution to the GZK [TJI2 puzzle. Due 
to their greater mass, their threshold energy for pion pho- 
toproduction is larger than for a proton. For this rea- 
son, an uhecron's energy degradation through the CMB 
is much smaller when compared to a proton jS] and it can 
come from farther away. A thorough search [9] for the 
source of the highest energy cosmic ray ever detected (by 
the Fly Eye's collaboration jTO]), pointed to a faraway 
[z = 0.545) quasar (3C147) as one of the best candi- 
dates. Although a proton coming from this distance can 
not reach the Earth, an uhecron can. 

Uhecron candidates arc found in extensions of the stan- 
dard model of particle physics. Heavier uhecrons (masses 
> 50 GeV) were excluded [11] as UHECR. Besides other 
reasons, the air showers they produce have their maxi- 
mum too deep in the atmosphere. Among the surviving 



candidates are the heavy gluino lightest supersymmetric 
particle (LSP) [T5J |T3j and strongly interacting wimpless 
dark matter [Hj. A search for the heavy gluino LSP us- 
ing CDF [13] and LEP data [15] constrained its mass to 
a 25 to 35 GeV window. Here we show that the neutral 
mode of this particle can be detected by UHECR tele- 
scopes, and this mass window allows for discrimination 
from the bulk of UHECR assuming it is composed by 
protons or nucleus. 

Our investigation [16 follows the uhecron definition 
stated in [11 . It is an electrically neutral, strongly in- 
teracting heavy exotic hadron. The bulk of its mass is 
carried by a single constituent. This is surrounded by 
hadronic degrees of freedom, which are responsible for 
the uhecron interaction. 

We simulate uhecron induced air showers in a simi- 
lar way as described in [IT] and then the detection and 
event reconstruction by a fluorescence detector (FD) as 
described in |17j . Proton and uhecron induced showers 
are compared and their discriminating parameters arc de- 
termined. As all UHECR simulations extrapolate known 
physics at lower energies to much larger energies, it is 
important to note that we compare uhecron to proton 
observables. In this way, we reduce the bias introduced 
due to uncertainties in the extrapolation of interaction 
models to high energies, since these uncertainties will af- 
fect both protons and uhecrons. 

As a result, we show that uhecrons with masses be- 
low 50 GeV can be detected in UHECR telescopes and 
discriminated against protons and nucleus. 

In the next section we describe our simulation of uhe- 
cron induced showers. Follows the description of the FD 
detection and event reconstruction simulation. We then 
describe the main uhecron induced shower features and 
compare them to proton and iron induced showers. In 



section IV we describe how to discriminate between pro- 
tons and uhecrons. In the following section we discuss 
photon induced showers. The last section presents our 
conclusions. 
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II. UHECRON INDUCED SHOWER 
SIMULATION 

When a UHECR impinges the Earth atmosphere, it 
generates a shower of particles. As the shower develops, 
the number of particles increases until it reaches a max- 
imum at a certain point in the atmosphere (Xmax). At 
this maximum the energy of each particle is low enough 
to be lost through ionization. The development of the 
shower as a function of the atmospheric depth (longitu- 
dinal profile) depends on the primary cosmic ray com- 
position. The longitudinal profile integrated energy is 
proportional to the primary cosmic ray energy. 

Air shower simulations include a particle cascade de- 
velopment integrated with an event generator. The later 
simulates the interactions between particles with air nu- 
cleus while the shower development simulates the parti- 
cle cascade versus atmospheric depth. In our simulation 
we use the Air Shower Extended Simulations (AIRES) 
package (v2.8.4a) [19] with SIBYLL (v2.1) [20JE] as the 
event generator. 

In order to simulate uhecron induced showers, we mod- 
ified both AIRES and SIBYLL. While modifications to 
AIRES are straight forward, and basically requires the 
inclusion of a new particle in the shower development, 
the modifications to the event generator are more com- 
plex. We use the modifications described in details in 

HQ. 

SIBYLL [20 models the interaction of a particle with 
an air nucleus as a combination of a low energy hadron- 
hadron interaction and a model for the "hard" part of the 
cross section. It also models hadron-nucleus interactions 
[21) . The interactions that occur high in the atmosphere 
have very large center of mass (CM) energies. SIBYLL 
extrapolates the known physics at much lower energies to 
higher CM energies (#(100 TeV)) using the dual parton 
model [24 superposed by minijet production. 

In short, the main modifications to SIBYLL (described 
in more details in |llj ) are as follows. The uhecron is rep- 
resented as a heavy single constituent (Q) surrounded by 
light hadronic degrees of freedom. Its interaction is sim- 
ulated in the same way as hadron-hadron interactions, 
which are represented |20j by production and fragmenta- 
tion of QCD strings. However, uhecron interactions use 
harder structure and fragmentation functions than the 
ones used for normal hadrons. In analogy to the B me- 
son, we describe the fraction of energy z carried by Q, 
using the Peterson fragmentation function [2"2"1 12"3"] : 



1 

z 1 



(2.1) 



where eq is proportional to Ag CD /mg. 

The good agreement between this fragmentation func- 
tion and data is described in (25]. This guarantees that 
most of the uhecron momentum is carried by the heavy 
constituent. The same function is used for the struc- 



ture function, which describe the fraction of energy of 
the hadron carried by Q. 

As the light constituents are responsible for the inter- 
actions, we take the uhecron-nucleon (<jjjn) cross section 
as the one for pion-nucleon interactions. Other modifi- 
cations are related to diffraction dissociation, where the 
lower mass limit of the excited state was modified accord- 
ing to the uhecron mass (mu). Also the "hard" part of 
the cross section with large momentum transfer, which 
is simulated as minijet production, is turned off for uhe- 
crons, since most of the momentum is carried by Q which 
does not interact. 

Figure [l] shows the average longitudinal profile of 320 
and 50 EeV iron, proton and uhecron (with mo = 20 
and 50 GeV) induced showers based on 500 showers for 
each primary. As uhecrons have less energy available for 
interactions than protons, its shower Xmax position is 
deeper in the atmosphere. As the uhecron mass increases, 
the available interaction energy decreases and the Xmax 
is deeper. These profiles show a fit with the Gaisser- 
Hillas [26 function (GH) to the simulated data. 
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FIG. 1: Average longitudinal profiles based on 500 iron, pro- 
ton and uhecron (with mu = 20 and 30 GeV) induced show- 
ers. Primary energies are equal to 320 EeV (top) and 50 EeV 
(bottom). These showers were generated at a 60° zenith an- 
gle. 

The Xmax position and number of particles at this 
position (Nmax) of each average profile in Figure [l] is 
shown in Table |T] The average Xmax position of 20 GeV 
uhecrons is about 100 g/cm 2 deeper than the one for pro- 
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TABLE I: Nmax and Xmax (slant depth) for shower profiles 
shown in Figure[T]and for a 30 GeV uhecron. Primary energies 
are 320 and 50 EeV. 



Energy (EeV) 


320 


50 


Particle 


Nmax 


Xmax 


Nmax 


Xmax 




(xlO 11 


) (g/cm 2 ) 


(xlO 10 ) 


(g/cm 2 ) 


Iron 


2.34 


797.1 


3.68 


749.4 


Proton 


2.23 


897.6 


3.58 


852.1 


Uhecron (20 GeV) 


1.94 


997.7 


3.06 


953.6 


Uhecron (30 GeV) 


1.92 


1005.3 


3.00 


967.4 


Uhecron (50 GeV) 


1.85 


1021.5 


2.90 


977.6 



tons for both primary energies. Although the longitudi- 
nal profile fluctuates, it already indicates that uhecrons 
resemble more protons than iron nucleus. For this rea- 
son we will determine ways to discriminate uhecron from 
protons. Our distributions show that the same procedure 
will more efficiently separate them from iron. 
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FIG. 2: Xmax distributions for 500 proton and 50 GeV uhe- 
cron induced showers of 320 (top) and 50 EeV (bottom). 

Although the average longitudinal profile gives an idea 
of the general differences between proton and uhecron in- 
duced showers, the individual profile fluctuates a lot. The 
Xmax distributions give a better idea of the fluctuation as 
well as the superposition among the different primaries. 
The Xmax distributions for proton and uhecron induced 
showers are shown in Figure [2j As uhecrons interact al- 
ways softly, the shower fluctuations are larger and the 
Xmax distribution is more spread than for protons. This 
feature will help in the proton uhecron discrimination. 



However, as we will show in the next section, showers 
with Xmax deeper than ground level are not accepted by 
the FD reconstruction. This requirement ends up lower- 
ing the uhecron acceptance. 



III. SIMULATION OF FLUORESCENCE 
DETECTION AND EVENT RECONSTRUCTION 

Fluorescence telescopes detect fluorescence photons 
emitted when charged particles transverse the atmo- 
sphere. As the air shower develops, the light emitted at 
different depths is collected by the FD photomultipliers 
(PMTs) and can be translated into an energy deposition 
longitudinal profile. The integration of this profile over 
the full shower path is proportional to the shower calori- 
metric energy. A small fraction (~ 10%) of the total 
shower energy is missed, since it is carried by neutrinos 
and by high energy muons which reach the ground. 

After generating uhecron, proton and iron induced 
showers, we simulate their detection by FDs. We use the 
same FD simulation as described in [17 , which followed 
the general procedure in [27]. As our simulation aims 
detection of rare events, a large coverage area is needed. 
For this reason we used the Pierre Auger FD parame- 
ters. The telescope altitude is set to 1500 m above sea 
level, 3.8 m 2 aperture covering an elevation angle from 
2° to 32° and using 1.5° pixel size PMTs. We take the 
telescope efficiency as 20%. 

In short, our simulation |17j translates the shower en- 
ergy deposited at each atmospheric depth into produc- 
tion of fluorescence photons. The propagation of these 
photons to the detector PMTs takes into account atten- 
uation due to Rayleigh (molecular) and Mie (aerosol) 
scattering [55] • Details of fluorescence detection such as 
effective collection area, mirror reflectivity, filter trans- 
mission, phototube quantum efficiency, noise and back- 
ground are included. Once the sequence of signals in each 
PMT is determined, we simulate the energy reconstruc- 
tion. We fold a 5° Gaussian error into the shower axis 
direction and transform back the PMT signal into de- 
posited energy. All effects that were taken into account 
in the detection simulation, are now determined by the 
new reconstruction shower geometry. 

We generated 2000 showers for each primary at 3 en- 
ergies (50, 100 and 320 EeV), all with a 60° zenith angle. 
Uhecrons with 20, 30 and 50 GeV mass were simulated. 
Each of these sets were used as inputs in the FD simula- 
tion. Each input was used 20 times, each with a different 
zenith angle and core position |17j . in order to simulate 
an isotropic flux. Overall, 40K FD events were simulated 
for each energy and particle. 

Once the longitudinal profile was reconstructed by the 
FD simulation a GH function was fit to determine the re- 
constructed energy. In order to cut badly reconstructed 
events, basic quality cuts were applied. These are listed 
in Table [IT] and are always apllied in our FD event recon- 
struction. All cuts but the GH fit x 2 are typically used 
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TABLE II: Quality requirements over FD simulated data. 
Events that do not meet these requirements are cut. All but 
the GH x 2 are found in 29 . "I> is the angle between the 
shower axis and the ground. 

Quality Requirements 



Hit PMTs 
Track length 
Zenith angle 
Xmax 
$ angle 
X 2 (GH fit) 



> 5 

> 200 g/cm 2 

< 60° 
visible 

< 132° 

< 50 



in Auger analysis [23]. The GH fit \ 2 was relaxed since 
this fit is not as good for uhecron longitudinal profiles 
as for proton's. $ is the angle between the shower axis 
and the ground and is used to minimize the Cherenkov 
contamination. 

Detection and energy reconstruction induces errors in 
the reconstructed longitudinal profile. Figure [3] compares 
the Xmax distributions for protons and 50 GeV uhecrons 
with 320 EeV primary energy, before and after the FD 
reconstruction. For a better visualization, we also show 
the distribution after FD reconstruction normalized to 
the number of input events. 
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FIG. 3: Xmax distributions before and after FD reconstruc- 
tion. Plots are for protons (top) and 50 GeV uhecrons (bot- 
tom) with 320 EeV primary energy. 

The large reduction in the number of events comes 
from geometrical factors as well as from the quality re- 



quirements. The shower detection is largely dependent 
on the inclination of the shower, core position and the 
detector field of view [TTJ [3D] . After the FD reconstruc- 
tion, both proton and uhecron distributions arc shifted to 
lower values and are also broader. The shift and reduc- 
tion of events is weaker for protons when compared to 
uhecrons. While detection uncertainties broadens both 
proton and uhecron distributions, the FD acceptance fa- 
vors lower Xmax values | 30l . For this reason more uhe- 
cron events are cut and the larger Xmax side of the dis- 
tribution is less accepted. This shifts the distributions to 
lower Xmax values. 

Since lower energy showers have Xmax at higher alti- 
tudes, they will be less affected by the FD acceptance. 
Our distributions follow this trend: for lower primary 
energies the Xmax distribution does not shift to lower 
values as much as for larger energies. Also, the reduction 
in the number of events is lower. While 84.7% (74.3%) 
of 320 EeV uhecrons (protons) are cut by the FD recon- 
struction, 81.5% (70.5%) of 100 EeV uhecrons (protons) 
are cut. 

Figure [4] shows the normalized maximum deposited en- 
ergy distributions (dE/dx)max before and after the FD 
reconstruction simulation for 320 EeV protons and uhe- 
crons (with 50 GeV mass). As for the Xmax distribu- 
tions, the maximum deposited energy also shifts to lower 
values. While the broading of the proton distribution due 
to detection and reconstruction errors is clear on both 
sides, the effect on uhecrons is not that clear, specially 
at lower deposited energies. This can be explained by the 
inherent uhecron shower characteristics, which fluctuates 
much more than proton showers. 

We also compare the reconstructed energy with the 
primary energy. The reconstructed energy is obtained 
by adding the missing energy to the calorimetric energy. 
While the latter is determined from the integration of 
the energy longitudinal profile, the missing energy is pa- 
rameterized from Monte Carlo simulations. We used the 
same missing energy parametrization as determined for 
protons in [3T] . Figure [5] shows the reconstructed energy 



error (given by (E rec - E P 



,)/E F 



7 ), before and 



after the FD reconstruction. An energy error of about 
3% is already observed in the reconstructed energy be- 
fore the FD simulation. This error is due to the missing 
energy parametrization, which was determined based on 
Corsika/QGSJET [32, 33 simulations and generates this 
error when using AIRES/SIB YLL. Our investigation will 
not be biased by this error, since we always compare uhe- 
crons with protons, and both are equally affected by the 
missing energy parameterization. 

As shown in the same Figure, the proton energy error 
peaks at the same energy before and after the FD simula- 
tion. While it is symmetrically distributed, the uhecron 
distribution is asymmetric. The main reason for this, is 
that the GH function is not the best fit for uhecron pro- 
files. Among other problems it does not account for the 
profile tail. In this analysis we did not attempt to find 
a better fit, but eventually it can help uhecron discrim- 
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FIG. 4: Normalized maximum deposited energy distributions 
before and after FD reconstruction. Plots are for 320 EeV 
primary energy protons (top) and 50 GeV uhecrons (bottom). 
The dashed curves represent the distribution after the FD 
simulation normalized to the number of events before the FD 
simulation. 



FIG. 5: Energy error distributions before and after FD recon- 
struction. Plots are for 320 EeV primary energy protons (top) 
and 50 GeV uhecrons (bottom). The green (proton) and blue 
(uhecron) narrow distributions are the errors before FD recon- 
struction, while the black (proton) and red (uhecron) broad 
distributions are the errors after the FD reconstruction. 



ination. As a result, uhecron showers will in average be 
reconstructed as lower energy showers, with a systematic 
energy error around -10%. 



IV. UHECRON - PROTON DISCRIMINATION 

As was shown in the previous section, the main char- 
acteristics of uhecron induced showers are larger Xmax 
with less particles at this position (lower Nmax) and 
slower development when compared to proton induced 
showers. Here we demonstrate the possibility of discrim- 
inating uhecron from proton induced showers using FD 
observables. Nucleus induced showers have even smaller 
Xmax and are more easily discriminated from uhecrons. 

Other than the Xmax and the (dE / dx)max (which has 
the same discriminating power as Nmax), the zenith an- 
gle 8 Z and the altitude Hmax at which the first light is 
detected by the FD can be used to discriminate uhecrons 
from protons. 

It can be seen from Figures [T] and [2] that uhecrons have 
deeper Xmax than protons. As a consequence a large 
fraction of uhecron induced showers that come vertically 
into the atmosphere are cut by the FD reconstruction. 
The requirement that the shower Xmax is visible (see 
Table fTTl) cuts most of the vertical uhecron showers. For 



this reason, uhecrons are better accepted at larger zenith 
angles and a cut on low 9 Z showers will be more effective 
on protons. 

As described in section [n] most of the uhecron energy 
is not available for interactions. For this reason, its first 
interaction point with a deposited energy larger than the 
FD threshold, will be deeper in the atmosphere than the 
first light collected from protons. Therefore Hmax can 
also be used as a discriminator. 

Figures [6] and [7] show distributions for the observables 
used as uhecron discriminators. All distributions are af- 
ter the FD reconstruction, for 320 EeV showers induced 
by protons and by 50 GeV uhecrons. 

In order to optimize all cuts on the discriminating pa- 
rameters (which from here on we call analysis cuts) , min- 
imizing the background contamination and maximizing 
the number of uhecrons, we use the following quality fac- 
tor: 



(N u + N p ) 



x JV» 



(4.2) 



where N u and N p are the number of uhecrons and protons 
after all analysis cuts were applied and a is a constant. 
Parameter a sets the strenght of the cuts. 

The quality factor q has to be maximized as a function 
of the analysis cuts. To achieve this maximization, we 
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FIG. 6: Xmax (top) and (dE/dx)m&x (bottom) distributions 
after FD reconstruction for 320 EeV primary energy proton 
and 50 GeV uhecron induced showers. The arrows show the 
position of the optimized analysis cuts. 



FIG. 7: 6 Z (top) and Hmax (bottom) distributions after 
FD reconstruction for 320 EeV primary energy proton and 
50 GeV uhecron induced showers. The arrows show the posi- 
tion of the optimized analysis cuts. 



scan q over a fixed range for each of the discriminating 
parameters Xmax, {dE/dx)vaax, Hmax and 6 Z . Each 
combination of cuts will yield a different q factor and 
the maximum value will indicate the optimized set of 
cuts. The arrows shown on Figures [6] and [7] indicate the 
cut values on each of these parameters. Parameter a in 
Equation 4.2| is set to 0.2. 

Figure 8] shows the Xmax distribution for 320 EeV 
showers generated by protons and by 50 GeV uhecrons 
before and after the analysis cuts were applied. It shows 
the effectiveness of the analysis cuts, since most of the 
protons are cut while a significant uhecron fraction sur- 
vives. The accepted region for each discriminating pa- 
rameter as well as the fraction of surviving events are 
shown in Table IIIII 

As discussed at the end of the previous section, uhe- 
crons will have their primary energy reconstructed with 
about a 10% error to lower values. For this reason we 
also compare 320 (100, 50) EeV proton showers with 
352 (108, 54) uhecron showers, corresponding to a 10% 
(8%) correction to the uhecron reconstructed energy. The 
results are shown on Tables |III| and |IV| where the first 
table compares 50 GeV and the latter 20 GeV uhecrons 
to protons. 

Figure [9] shows both Xmax and (dE/dx)max distribu- 
tions for 320 EeV protons and 50 GeV uhecrons with 
320 EeV and 352 EeV primary energy. The Xmax 
distribution will not change significantly although it 
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FIG. 8: Xmax distribution for protons and 50 GeV uhecrons 
of 320 EeV primary energy before and after analysis cuts were 
applied. Cut values and fraction of events surviving the cuts 
are shown in Table UTTl 



shifts slightly to deeper Xmax values. However the 
(dE/dx)max. distribution changes significantly since the 
energy correction implies in a larger deposited energy. 
This shift in the (dE/dx)max distribution will reduce 
the discriminating power of this observable. 

Both uhecron energy and mass are important factors 
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11.50 


973.3 


108 


100 


0.366 


0.039 


0.143 


157e+08 


0.637 


11.44 


956.3 


54 


50 


0.299 


0.016 


0.080 


6.64e+7 


0.400 


11.41 


882.8 



TABLE III: Fraction of events after analysis cuts. E v and E u are primary energy (in EeV) of protons and 50 GeV uhecrons, 
respectively; N p (N p o), N u (N u o) and Nt are the number of protons; uhecrons and the sum of proton with uhecron induced 
showers after the FD simulation and after (before) all analysis cut are applied. The last 4 columns indicate the accepted region 
for each discriminating parameter after cut optimization, in units of GeV cm 2 /g; rad; km and g/cm 2 , respectively. 
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TABLE IV: Same as Table 11111 but now uhecrons have 20 GeV mass. 
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FIG. 9: Xmax (top) and (dE/dx)ma,x (bottom) distributions 
after all analysis cuts were applied for 320 EeV protons and 
50 GeV uhecrons with 320 EeV and 352 EeV primary energy. 



for discrimination from other primary particles. An en- 
ergy increase favors deeper Xmax parameters and there- 
fore less FD acceptance but on the other hand enhances 
the intrinsic differentiating shower characteristics in re- 
lation to protons. Lower uhecron masses approximate 
their shower intrinsic characteristics to proton showers, 
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0.05 


0.6 
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0.0015 
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54 
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0.01 


0.8 


0.081 


0.0001 


0.156 



TABLE V: Fraction of events after analysis cuts for protons 
and 50 GeV uhecrons. <j> u /4>p is input ratio of uhecron to pro- 
ton induced showers into the FD simulation, a is a parameter 
in Equation |4.2| and all other parameters are described in Ta- 
ble EI] 



but increases the uhecron FD acceptance. As shown in 



Tables III and |IV| it is possible to greatly reduce the pro- 
ton contamination using Xmax, (dE / dx)max, Hmax and 
Z as discriminating parameters. The proton contamina- 
tion in the final event sample is at maximum 15% for 
50 GeV uhecrons and below 20% for 20 GeV uhecrons. 



A. Uhecron — proton flux ratio 

Up to now we have considered the same number of 
input uhecron and proton induced showers into the FD 
simulation. This is equivalent to an equal flux of protons 
and uhecrons arriving at the Earth. However, consider- 
ing the latest UHECR flux measurements [3J H] a much 
lower uhecron flux has to be considered at least up to 
energies around the expected GZK cutoff. Beyond this 
point, a nucleon or nucleus flux is not expected. Events 
at energies beyond the GZK cutoff might indicate new 
physics. 
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Here we redo the analysis described in section IV but 
reducing the uhecron flux 
to the proton flux 



to 10, 5 and 1% relative 
r p . We analyse a 1% uhecron frac- 
tion for a 50 EeV primary energy (54 EeV for a uhecron 
shower) since at lower energies uhecrons might be present 
as a small fraction of the flux. Even in this scenario it is 
possible to discriminate uhecrons from protons. We sum- 
marize our results in Table [Vj In order to enhance the 
final number of uhecron events, we use larger a parame- 
ter values (see Equation 4.2 1. After all analysis cuts are 



applied, the proton contamination in the final sample, 
for a 1% uhecron flux is around 16%. As we will discuss 
in the last section, our results indicate the feasibility of 
discriminating uhecrons from protons, even with a much 
smaller uhecron flux. 



B. Sample independence test 

In order to check our uhecron analysis and the dis- 
criminating power of the Xmax, (dE / dx)max, #max and 
Hmax observables, we applied the same analysis cuts as 
described above to a new set of data. This new set of data 
uses the same 2000 showers generated from our shower 
simulation, for each primary particle (where 3 different 
uhecron masses - 20, 30 and 50 GeV - were assumed) 
and for each different primary energy (50, 100 and 320 
EeV) and input it with different geometry [T7] than the 
original analysis to the FD simulation. The analysis cuts 
applied to this new simulated data set were the ones de- 
termined in the original analysis. 

We obtain similar results as in the original analysis. 
The analysis cuts have the same discriminating power. 



Tables |VTJ and [VII] summarize the analysis results for this 
new data set (for 50 GeV and 20 GeV uhecrons respec- 
tively). As can be seen the results are compatible with 
the ones in Tables UTTI and Hvl 



V. UHECRON - PHOTON COMPARISON 

Photon showers develop deeper in the atmosphere than 
proton showers. For this reason it resembles more uhe- 
cron than proton induced showers. However it is impor- 
tant to note that the competition among uhecron and 
photons is not realistic. At these energies, both photons 
and uhecrons are proposed in beyond the standard model 
of particle physic scenarios. These either propose ultra 
high energy photons or exotic hadronic particles. It is 
already known that the photon fraction of the UHECR 
flux is very small, which constrain many top-down mod- 
els [34] . In these models, ultra high energy photons would 
be produced from exotic heavy particle decay. Auger re- 
sults limit [IB] the photon fraction of the UHECR flux to 
2% (5.1% and 31%) of the total flux above lx (2x and 
4x) 10 19 eV with 95% CL. 

It is also important to note that photon induced show- 
ers develop differently from hadronic induced showers. 



This difference is maily due to smaller particle multiplic- 
ity in the eletromagnetic cascade when compared to a 
hadronic cascade. As a consequence the photon shower 
Xmax is in average deeper than the proton Xmax. Also, 
the number of muons in hadronic showers is greater than 
in photon showers, due to charged pion decay. For this 
reason, a FD uhecron - photon discrimination can be 
greatly enhanced by ground detector information. 
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FIG. 10: Longitudinal profile for photon induced showers 
with and without the preshower effect. Proton and 50 GeV 
uhecrons are also shown. Proton and photon showers have 
320 EeV primary energies while uhecron has 352 EeV. All 
profiles are before FD reconstruction. 

Another important effect to be taken into account is 
due to photon interactions with the Earth geomagnetic 
fields. As a consequence they preshower before entering 
into the atmosphere. In Figure 10 we compare proton, 



uhecron (with 50 GeV mass) and photon longitudinal 
profiles for 320 EeV induced showers. Photon profiles are 
shown with and without preshower. The effect seen in the 
longitudinal profile will be present in the Xmax distribu- 
tion as well. The preshower effect changes both photon 
longitudinal and Xmax distributions in a way that it will 
resemble more uhecron showers than if this effect was 
not present. However as mentioned above the hadronic 
characteristics of uhecron showers might allow for their 
separation [33] , 



VI. DISCUSSION AND CONCLUSION 

We have shown that UHECR experiments, such as the 
Pierre Auger Observatory, have the potential to detect 
exotic massive hadrons. Also that it is possible to dis- 
criminate them from nucleon induced showers. 

Although both proton and uhecrons produce hadronic 
showers, uhecron characteristics will allow discrimination 
from protons. As the uhecron mass increases, its induced 
shower develop more slowly and fluctuates more. These 
characteristics allow to better distinguish heavier uhe- 
cron from proton showers, although it is also possible 
to distinguish lighter uhecrons. While the proton con- 
tamination in the final simulated data sample, after all 
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TABLE VI: Same as Table [III] Cuts are now applied to new data set. Accepted region for each discriminating parameter was 
defined from original data set. Uhecron mass was set to 50 GeV. Results with original and with new data set are compatible. 



E u 


E p 


N u /N u0 


Np/N p o 


N p /Nt 


(dE / dx)ma,x 


0z 


Hmax 


Xmax 












> 


> 


< 


> 


352 


320 


0.390 


0.062 


0.198 


5.54e+08 


0.712 


11.41 


961.4 


108 


100 


0.359 


0.057 


0.188 


1.74e+08 


0.616 


10.85 


951.7 


54 


50 


0.411 


0.071 


0.198 


8.12e+07 


0.300 


10.90 


922.3 



TABLE VII: Same as Table ED but now uhecrons have 20 GeV mass. 



analysis cuts are applied, is at maximum 15% for 50 GeV 
uhecrons it is around 20% for 20 GeV uhecrons. 

We have also shown the effects of fluorescence detection 
and event reconstruction. FD requirements can exclude 
uhecron showers that are naturally better discriminated 
from protons. However even after FD detection and event 
reconstruction it is possible to separate showers induced 
by these two primaries. We have shown that FD observ- 
ables such as Xmax, (dE / dx)ma,x, #max and Hmax, are 
good discriminators. 

Although we have no prediction for the ratio between 
proton and uhecron induced showers, we have shown that 
the uhecron flux can be as small as 1% of the total flux 
and still be discriminated from protons. At lower primary 
energies, standard model particles should dominate the 
UHECR spectrum, whereas at energies beyond the GZK 
cutoff it is possible to have a larger exotic flux. 

It is important to note that the search for beyond 
standard model particles is complementary to acceler- 
ator searches. It depends on an assumed prior model. 
If the Large Hadron Collider (LHC) has indication of a 
heavy gluino [TSJ [T3] , UHECR telescopes can look for it 
in a complementary way. Or vice- versa, one can find uhe- 
cron candidates among UHECR and depending on LHC 
results investigate its identity. Heavy gluino |12[ 113] and 
strongly interacting Wimpless particles [2] are examples 
of uhecrons. The current allowed heavy gluino mass win- 
dow [T3] (25 to 35 GeV) is within the uhecron mass limits 



that allow separation from proton or nuclei background. 

We have also shown that our method for uhecron detec- 
tion and background reduction is independent from our 
simulated data. After our analysis method was deter- 
mined, we applied it to a new data set. The new results 
show that our discriminating parameters have the same 
power as when applied to the original data set. 

We also have compared uhecron to ultra high energy 
photon induced showers. Although it is not expected to 
have both these particles as UHECR primaries, the dif- 
ferences between a hadronic and photon induced shower 
should allow for their discrimination [35] ■ Ground detec- 
tors should improve this discrimination. 

As the uhecron flux at energies around the GZK cut- 
off has no reason to be large, the construction of the 
northern Auger site will dcfinctcly improve the uhecron 
detection probability. 
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